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We propose to determine the coefficients in the Fermilab heavy quark action by matching the finite-volume, 
off-shell, gauge-fixed propagator and vertex functions with those determined in the exact, relativistic theory. 
The matching relativistic amplitudes may be determined either from short-distance perturbation theory or from 
finite- volume step-scaling recursion. Tree-level matching and non-perturbative off-shell amplitudes are presented. 



1. Introduction 

The Fermilab action PP is the least restric- 
tive framework permitting lattice methods to be 
used for heavy quarks when the quark mass is 
large compared to the inverse lattice spacing, 
TOquarka > 1- Howcver, at present the coefficients 
of this action have only been calculated using 
low order lattice perturbation theory. In order to 
have better control of errors for the physical quan- 
tities calculated on the lattice, we propose to de- 
termine these coefficients non-perturbatively. We 
plan to achieve this by imposing a matching con- 
dition between the off-shell quark propagator and 
quark-gluon vertex function determined from the 
Fermilab action, and the same propagator and 
vertex computed under more physical conditions. 

We propose two ways to achieve this goal. The 
most naive method would be matching the prop- 
agator and vertex function computed in the con- 
tinuum theory up to some order in perturbation 
theory with the RI-MOM NPR calculation 
on the lattice (Figure . This would be done 
through the "matching window" I^: Ag^^ <C 
— ^ 1/a? . The alterative will be purely 
non-perturbative, by connecting the original cal- 
culation with one performed at a smaller lat- 
tice spacing a' = a * e (as shown in Figure EJ . 
A further advantage of such an entirely non- 
perturbative matching is that the kinematic "win- 
dow" for performing such matching now has no 
lower bound since there is no requirement that 
the perturbation theory be accurate. 



*We thank RIKEN, Brookhaven National Laboratory and 
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Lattice Continuum 

Figure 1. Proposal-I. Match the lattice result to cor- 
responding continuum perturbation theory calcula- 
tion to sufficient order, in the usual RI-MOM scheme. 
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Figure 2. Proposal-II. Match the lattice re- 
sult, through a step-scaling technique, to finer non- 
perturbative, 0(a) improved light quark calculations. 

2. Off-shell Fermilab action 

Extending the off-shell 0(a) improvement of 
light quarks 0] to heavy quark, we will adopt the 
approach of using an unchanged fermion action 
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in the simulation: 

5* = ^ IpniDpL + m)tpn 

and improved quark fields of the form: 
iP^iP^ Zg^'[l + aCq'iDpL + m) 

- - — i - 

4i ^ Zq + acg'(-£'jr^ + m) 

where Dp^ is the Fermilab Dirac operator. Note: 
Zq and c's are all functions of ma. In summary, 
we have eight parameters as a function of ma to 
be determined: fg, ce, cb (in the action) and 
Zq, Cq , c%Qj and c%qj (in the improved quark 
fields). 

If we look at the low spatial momentum behav- 
ior of the Dirac structure we find that there are 
more than the six required conditions from 
and the two from S{p) , and thus we are able to 
fix these eight parameters. 

3. Tree-level coefficient calculations 

Before we start the complicated computer cal- 
culations, we can test our proposal-I by starting 
with lattice perturbation theory in tree-level and 
matching to the corresponding continuum quan- 
tities: 

= Zq{ii; ay^Z^^^{n; a)Ajf*(p, q; a)\„^2_p2^^2 
{{S{p))-'U2_p2^^,2r'^y'^ 

= tZq[ipoj,A'f{p'^) + ^p.n^A'r{p')] + B'^'ip') 
Z3D{q'^;fi)\q2^ 2 = — 

Expanding in the variables pia and q^a up to 
0(a) we determine: 

C = cos(poa + 9oa/2) = cos(po«) — goa/S sin(poa) 
Z^^ = (1 + 4:Cqma)cos{poa + qoa/2) 
t _ 1 r_Poa_ 
2ma smpoQ. 

X (cos(poa + 90^1/2) — 1)(1 + 4cgma)} 

Cq = 1/2 - c%(ji cos(poa + (7oa/2) 
ce = 2cg 
Cb = 2c, C 



Taking the p^a and ma ^ 1 limit, we reproduce 

the known coefficients 

CE = CB = C = rs = 1, c^gI = , 

Cq = 1/2, Zq ~ 1 — 2ma 

for the improved, light quark action. 

4. Non-perturbative Green's functions 

We have studied connected 2-point Green's 
functions (quark and gluon propagators) 6 and 
have produced answers consistent with published 
papers. The accuracy achievable with moderate 
statistics is adequate to determine those improve- 
ment coefficients that are constrained by the 2- 
point functions. From now on, we will focus on 
the least known quantity: the quark-gluon vertex 
function. 

The quark-gluon vertex function can be calcu- 
lated from the following 3-point function j^]: 

The Fourier transformed amputated vertex func- 
tion in momentum space is 

= {S{p))~'V:{p,q)t,,{S{p + q))-'{D{q),,)-' 

We calculate the quantity jIm T'''(7p ^JTa) in 
the case where q = 0. 

As a first step we plan to non-perturbatively 
determine the coefficients in light quark clover 
action, since only four parameters are involved: 
Csw (in the action) and Zq, Cq , and CNGiij^ the 
improved quark fields). Our gauge configurations 
use the Wilson gauge action with j3 = 6.0 , a~^ = 
1.922 GeV, 16^x32 lattice size. 

We use Landau gauge and proceed by two ap- 
proaches: 1. fix to maximal axial gauge (MAG) 
before Landau gauge fixing; 2. fix only to Landau 
gauge. We found that the two different proce- 
dures agree for the gluon propagator, have slight 
effects on quark propagators (decreasing as the 
momentum increases) but give very different re- 
sults in the vertex calculation [Hj, as shown in 
Figure 13 We believe the differences are caused 
by Gribov copies. From now on, we will follow 
procedure 2. for better control over which gauge 
orbit is selected by the Landau gauge fixing pro- 
cedure, and for the cleaner signal. 
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Figure 3. Gluon "q = 0" vertex function calculation 
from the tadpole improved clover action with k — 
0.137 from 43 configurations. 

Figure 0] shows the comparison of the unrcnor- 
mahzed q = vertex functions from NP clover 
with K = 0.13445 without off-shcU improvement, 
and with domain waU fermions (L5 = 16, M5 
= 1.8). DWF is automatically 0(a) off-shell im- 
proved and the vertex nicely becomes indepen- 
dent of momenta in high momentum limit, as one 
would expect from a perturbative analysis. It is 
precisely this difference in behavior that will en- 
able us to determine the coefficients to off-shell 
improve the clover action, first in the light quark 
limit, and later as a function of ma. The masses 
of the pseudoscalar meson calculated from these 
two actions are similar, around 500 MeV. Once 
we impose the 0(a) off-shell improved quark field 
clover action, the difference should be reduced. 

The "q = 0" kinematics has cleaner signal than 
"q ^ 0" and would be a better candidate for a 
step-scaling proposal, while the "q ^ 0" would 
be required to implement proposal-I 0. This 
would require better control of the statistical er- 
rors. The large error bar at "q 0" comes mainly 
from the weak correlations between the quark 
propagator and the gluon fields. We are currently 
working on using sources with definite momentum 
to achieve improved statistical accuracy. 

5. Summary and Outlook 

In this paper, we have proposed two meth- 
ods for non-perturbatively determining the co- 
efficients of the heavy quark action. We calcu- 
lated the tree-level coefficients and reproduced 
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Figure 4. Gluon vertex function with "q = 0" kine- 
matics for 379 configurations. 

the known light quark limit. It seems possible 
to determine the Fermilab coefficients with suffi- 
cient effort through either of our proposals. We 
will next look at other Dirac projections and mo- 
mentum configurations. Once we have investi- 
gated the optimal method more thoroughly, we 
will start with non-perturbative determination of 
Csw for the light quark action and reproduce the 
Alpha collaboration result. We will then move on 
to the large quark mass case: the Fermilab action. 
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